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A  study  of  the  effects  on  the  population  failure  rate  of  component 
failures  during  the  period  of  transition  to  steady  state  is  presented. 
Comparisons  between  the  simulated  failure  rates  and  the  system  failure 
rate  predicted  under  the  present  Air  Force  assumption  of  an  underlying 
Poisson  process  are  accomplished  for  various  MTBF's  and  parameters 
of  the  Weibull  distribution.  The  decreasing  hazard,  constant  hazard, 
and  increasing  hazard  functions  are  studied  and  presented. 

This  study  was  performed  as  a  simulation  of  the  transition 
period  using  the  GASP  II  package  developed  by  Pritsker  and  Kiviat. 

Such  a  simulation  package  greatly  simplifies  the  required  programs 
necessary  to  be  written.  The  compvler  program  used  in  this  study  is 
included  as  part  of  this  paper. 

The  results  of  the  study  indicate  there  are  significant  differ¬ 
ences  between  the  pi  esent  Air  Force  methods  of  predicting  failure 
rates  and  the  simulated  failure  rate  distribution  during  the  transition 
period.  Plots  of  the  simulated  results  compared  to  the  present  Air 
Force  methods  indicate  these  differences  co  ild  reach  1 0 C';l>  of  the  Air 
Force  estimated  failure  rate  when  the  Weibull  shape  parameter  is 
greater  than  3  and  the  components  j\i"r”F  exceeds  1000  hours. 

The  study  indicates  that  component  failure  rates,  which  can  be 
described  bv  a  Weibull  density  function,  will  result  in  oscillations 
about  the  present  Air  Force  computed  rate  at  steady  slate. 
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CHAPTER  I 


INTRODUCTION 


Structure 

For  the  convenience  of  the  reader,  we  have  structured  the 
following  report  so  that  Chapters  III  and  IV  are  not  required  for  an 
intuitive  understanding  of  the  background  and  results  of  our  investiga¬ 
tion.  A  reader  with  little  time  or  desire  for  mathematical  background 
material  may  wish  to  omit  Chapters  HI  and  IV  in  the  interest  of  brevity. 
In  addition,  we  offer  our  definitions  of  the  following  commonly  used 
terms. 

Component-  -The  general  term  for  any  part,  subassembly,  line 
replaceable  unit,  system,  or  piece  of  equipment  whose  failui’e  and 
replacement  must  be  considered  in  logistics  plans. 

Component  Population- -The  total  number  of  a  component  type  in 
the  Air  Force  operating  inventory. 

Failure--The  opposite  of  "operating  normally.  "  Thus,  for  our 
purposes,  a  component  has  only  two  states  of  nature;  either  it  is 


serviceable,  or  it  has  failed. 
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Failure  Rate- -The  rate  of  failure  in  a  population  of  identical 
components.  This  measures  applies  when  components  are  repaired  (to 
new  condition)  or  replaced  (by  new  components)  and  thus  continue  to 
operate.  It  is  computed  by  dividing  the  total  number  of  observed 
failures  in  a  time  period  by  the  accumulated  operating  time  during  the 
period. 

Mean  Time  Between  Failures  (MTBF)--The  average  operating 
time  expected  between  failures  in  a  population  of  identical  components. 
This  measure  has  meaning  only  when  we  are  discussing  a  population 
where  there  is  repair  or  replacement. 

Hazard  Rate- -The  instantaneous  failure  rate  of  the  components 
under  consideration;  or,  the  expected  failure  rate  during  the  next  time 
interval. 

Mean  Time  To  Failure  (MTTF)--The  average  operating  time 
expected  before  failure  of  a  component  which  is  not  repaired  or  replaced. 
This  is  simply  the  average  time  to  failure  of  "n"  units,  i.  e.  ,  the  sum  of 
"n"  individual  unit  times  to  failure  divided  by  "n"  units. 

Statement  of  the  Problem 

One  of  the  basic  factors  associated  with  the  logistical  support  of 
a  weapons  system,  or  a  group  of  weapons  systems,  is  the  determination 
of  future  requirements.  Specifically,  we  need  to  know  how  often  parts 
and/or  equipment  will  fail.  In  practice,  the  occurrence  of  a  component 
failure  is  so  basic  to  the  support  function  that  its  importance  is  so;  le- 
times  ovt  rlooked.  Yet,  the  frequency  and  distribution  of  component 
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failures  determine  how  many  spare  parts  are  to  be  purchased,  how 
many  support  manhours  are  to  be  scheduled,  what  base  and  depot 
support  facilities  are  to  be  made  available,  and  even  how  many  aircraft 
are  to  be  required  to  perform  the  mission. 

In  the  established  procedures  for  determining  support  require¬ 
ments  (5),  the  Air  Force  makes  an  implicit  assumption  as  to  the  under¬ 
lying  distribution  of  failures.  We  assume  each  component  population 
has  a  constant  failure  rate  with  respect  to  time.  Under  such  an  assump¬ 
tion,  we  ignore  the  reality  of  the  situation  where  components  actually 
deteriorate  over  time  and  that  the  probability  of  failure  may  increase  as 
the  components  accrue  larger  and  larger  operating  times.  Thus,  under 
such  an  assumption  of  a  constant  failure  rate,  we  may  unintentionally 
introduce  serious  errors  into  our  logistics  planning  to  meet  future  sup¬ 
port  requirements. 

Objective 

The  objective  of  this  research  is  to  investigate,  through  digital 
computer  simulation,  the  distribution  of  failure  rates  during  the  transi¬ 
tion  period  from  the  initial  acquisition  of  new  weapons  systems  to  the 
attainment  of  a  'steady  state"  condition,  resulting  from  the  incremental 
introduction  of  a  fleet  of  aircraft.  For  simplification  we  assume  each 
aircraft  consists  of  a  self-contained,  independent  component  whose 
probability  of  failure  can  be  approximated  by  a  Weibull  distribution. 
Specifically,  we  were  interested  in  the  change  of  the  component 
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population's  failure  rate  with  respect  to  time.  7n  reliability  terms,  we 
wanted  to  investigate  the  pooled  output  of  N  superpositioned  renewal 
processes  where  N(t)  =  £(t).  ^ 

Hypothesis 

The  fleet  wide  distribution  of  failures  for  any  replaceable  equip¬ 
ment/component  is  at  all  times  a  Poisson  process;  therefore,  the 
assumption  of  a  constant  failure  rat®  is  always  valid  in  estimating 
future  expected  failures. 


^We  are  interested  in  the  system  average  failure  rate  at  a  point 
in  time.  This  failure  rate  is  the  sum  of  the  individual  failure  rates 
divided  by  the  number  of  failures,  i.  e.  , 

n 

L 

N(t)  =  1=1  =  f(t) 


n 


CHAPTER  H 


BACKGROUND 


The  Air  Force  Method  for  Estimating  Part  Failures 

If  each  part,  subassembly,  or  line  replaceable  unit  (LRU)  in 
cur  supply  inventory  came  with  a  "service  life"  stamped  next  to  the 
part  number,  logistics  plans  for  future  requirements  could  be  made 
by  anyone  with  a  desk  calculator.  Components  could  be  replaced 
before  failure,  stocked  exactly  in  the  right  amounts,  and  always 
reordered  in  the  most  economical  quantities.  Unfortunately,  this  is 
not  the  case,  since  individual  component  times  to  failure  occur  randomly 
and  are  not  exactly  predictable.  To  overcome  this  problem  and  prepare 
meaningful  quantitative  support  plans,  the  Air  Force  has  adopted  a 
method  for  estimating  the  service  time  of  components  by  averaging 
service  times.  This  paper  is  concerned  with  the  accuracy  of  such  an 
estimation  method  during  the  initial  provisionin'?  process,  and  the  early 
life  of  a  weapon  system. 

The  established  Air  Force  procedures  for  predicting  component 
failure?  are  so  basic  that  we  seldom  realize  we  are  using  a  "method,  " 
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and  that  there  are  requirements  which  must  be  met  before  applying  it. 
The  following  is  a  rough  description  of  present  Air  Force  procedures 
(for  mathematical  derivations  and  references  see  Chapter  III). 

In  order  to  sensibly  deal  with  random  failures,  we  first  must 
make  a  few  assumptions.  One  o[  the  most  helpful  assumptions  is  to 
assume  that  while  Individual  components  fail  at  random,  the  total 
parent  population  is  in  a  steady  stale  condition.  For  our  purposes,  a 
component's  supply  system  would  be  in  steady  state  if  the  average  Air 
Force -wide  demand  for  the  component  remained  constant  year  after 
year.  Where  the  parent  population  is  in  steady  state,  or  sufficiently 
close  so  that  we  could  assume  steady  state,  then  we  can  use  the  popu¬ 
lation's  long  run  average  demand  as  an  estimator  for  a  future  period  of 
time.  For  example,  we  might  predict  1125  F-4  tire  failures  next  year 
because  we  have  averaged  1125  failures  for  the  last  five  years.  Such 
a  quantitative  measure  belongs  to  a  parent  population  in  steady  state 
and  identifies  the  average  number  of  failures  experienced  during  the 
prior  time  period. 

Next,  we  assume  that  component  failures  are  distributed  uni  ¬ 
formly  over  the  entire  interval.  If  we  have  a  homogeneous  system 
which  is  not  subject  to  large  seasonal  demands,  then  this  assumption 

^We  would  not  expect  a  particular  F-4E  aircraft  to  use  exactly 
the  same  number  of  right  main  tires  every  year;  however,  if  the  yearly 
F-4  flying  program  did  not  change,  we  could  expect  the  fleetwide 
demand  for  F-4  tires  to  be  approximately  cons‘ant. 
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will  be  valid.  We  can  now  calculate  an  estimate  of  the  population 

failure  rate  for  any  desired  small  time  interval.  Thus,  if  we  averaged 

1125  F-4  tire  failures  each  year  while  logging  225,000  flying  hours, 

we  cart  assume  we  experience  .  005  failures  per  flying  hour,  or  an 
« 

average  of  200  flying  hours  between  failures.  Since  the  component 
population  is  in  steady  state,  is  homogeneous,  and  composed  of  identi¬ 
cal  items,  we  can  assume  the  above  estimate  applies  uniformly  to  the 
individual  components.  Thus,  each  component's  failure  rate  and  time 
between  failures  would  be  identical  to  the  population  failure  rate  and 
time  between  failures,  in  our  example,  we  would  assume  each  F-4 
tire  had  a  failure  rate  of  .  005  failures  per  flying  hour,  or  an  average 
service  life  (MTBF)  of  200  flying  hours. 

The  procedure  is  complete;  a  quantitative  measure  has  been 
determined  which  can  be  used  in  estimating  future  requirements.  How¬ 
ever,  this  method  provides  quantitative  estimates  only  for  the  steady 
state  conditions.  Thus,  future  requirements  may  be  estimated  from 
current  supply  demand  rates.  As  the  component  and  its  parent  popu¬ 
lation  age,  these  estimates  can  be  refined  as  more  data  is  collected. 

Identical  or  Similar  Components 

When  we  begin  the  planning  process  for  a  new  weapons  system, 

<1 

■'Seasonal  demands  for  items  by  their  very  nature  do  not  exhibit 
a  uniform  distribution  throughout  the  time  interval.  However,  where 
the  seasonal  demands  are  relatively  small  in  comparison  to  total 
demands,  then  such  demands  would  not  significantly  affect  this  assump¬ 
tion. 
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we  do  not  have  to  start  over  for  every  new  component.  Many  items  on 
a  new  aircraft  are  identical  or  quite  similar  to  components  used  before. 
If  we  have  collected  data  or.  similar  systems,  we  may  be  able  to 
accurately  predict  future  requirements  and  provision  for  the  new  com¬ 
ponents  using  our  past  demand  data.  However,  if  the  basic  assump¬ 
tions  we  made  while  applying  the  previous  "method"  are  not  valid  for 
the  new  component  population,  sizable  errors  may  be  introduced  into 
provisioning  plans. 

If  the  new  aircraft  are  to  be  acquired  in  large  numbers,  say, 

100,  500,  or  700  in  the  total  fleet,  it  is  unlikely  that  all  the  aircraft 
will  be  placed  into  operation  immediately.  They  will  probab'v  be  intro¬ 
duced  into  the  active  inventory  over  a  period  of  months,  or  even  years. 
Does  our  first  assumption  of  a  steady  state  system  still  hold?  If  the 
components  we  are  buying  are  not  yet  members  of  a  steady  state  popu¬ 
lation,  can  we  attribute  to  each  of  them  the  old  system  failure  rate  and 
mean  time  between  failures?  In  some  cases  yes;  in  some  cases  no. 

To  illustrate  this  point,  we  will  use  a  hypothetical  situation. 
Consider  a  "specialist"  shop  in  a  maintenance  squadron  that  supports 
ten  aircraft.  This  shop  repairs  a  photographic  system  that  has  two 
components:  a  power  supply  unit  and  a  camera  unit.  Assume  that  the 
system  has  been  in  the  Air  Force  inventory  for  a  few  years  and  that  it 
is  in  steady  state.  The  specialist  repairman  is  new  and  eager.  On  his 
first  day  in  the  shop  he  erects  a  four-month  wall  chart  for  each  com¬ 
ponent,  complete  with  aircraft  tail  numbers  and  the  date.  For  four 
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months  he  religiously  records  each  failure  as  a  large  "X"  on  'he  chart 
until,  at  the  end  of  April,  the  charts  appear  as  in  Figure  1.  Tha 
charts  appear  quice  similar  with  five  failures  each  month.  Knowing 
that  each  aircraft  flies  100  hours  a  month,  our  camera  specialist 
calculates  each  component's  failure  rate 


5  failures 
1000  flying  hours 


.  005  failures  per  hour 


and  mean  time  between  failures 


1000  hours  ,AA  , 

-  =  200  hours. 

5  failures 

While  the  results  from  the  charts  appeared  equal  for  the  power 
supply  units  and  the  cameras,  they  disagreed  with  the  repairman's 
intuition.  He  knew  the  components  were  entirely  different;  and  there¬ 
fore,  he  guessed  that  they  should  possess  different  failure  patterns. 

The  power  supply  unit  was  a  solid  state,  fully  electronic  assembly  that 
appeared  to  fail  at  random.  Conversely,  the  camera  unit  was  com¬ 
posed  of  small  gears,  bearings,  and  sprockets  which  became  clogged 
with  dust-like  film  emulsion  and  thus  failed  at  regular  intervals.  To 
test  his  conjecture,  the  repairman  went  back  to  his  charts  and  con¬ 
nected  the  •  X"s  for  each  aircraft  with  a  solid  line.  The  lines  represent 
the  operating  time  prior  to  failure  for  each  individual  component 
(Figure  2). 

For  the  power  supply  components,  the  service  time  lines  were 
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many  different  lengths;  but,  the  average  length  was,  as  computed,  200 
hours.  The  probability  that  they  will  fail  in  the  next  moment  is  not  a 
function  of  operating  age.  Because  of  this  property  the  length  of  serv¬ 
ice  varies  randomly  from  very  short  to  very  long.  The  service  times 
for  the  camera  units  appear  very  consistent,  varying  only  over  a  small 
range  of  times  around  200  hours. 

Now  that  we  have  accumulated  failure  data  in  our  hypothetical 
shop,  suppose  v'e  go  one  step  further  and  attempt  to  provision  for  an 
entirely  new  camera  system.  Assume  the  squadron  acquires  ten  air¬ 
craft  from  the  production  line  equipped  with  the  new  camera  system 
and  expects  to  fly  each  aircraft  an  average  of  100  hours  a  month 
beginning  1  January.  We  use  our  hypothetical  (steady  state)  failure 
rates  and  stock  parts/equipment  based  upon  the  prediction  of  five  power 
supply  and  five  camera  failures  each  month.  Figure  3  shows  the 
hypothetical  results  experienced  on  the  basis  of  our  earlier  assump¬ 
tions. 

The  power  components  appear  in  steady  state  immediately,  how¬ 
ever,  the  camera  components  experience  no  failures  for  the  first  two 
months.  Then  they  begin  a  two  month  failure  cycle  which  continues 
until  the  small  variations  in  service  time  smooth  the  distribution  to  the 
steady  state  condition  shown  in  Figure  1.  We  have  over  stocked  parts/ 
equipment  for  the  camera  unit  during  the  months  of  January  and 
February  since  no  failures  were  experienced.  However,  during  the 
month  of  March,  we  discover  that  we  do  not  have  enough  parts/ 
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equipment  on  hand  to  satisfy  all  failures  and  the  mission  is  affected  due 
to  lack  of  parts.  The  stockage  requirements  for  repair  support  of  the 
camera  system  may  be  increased  to  insure  sufficient  spares  are  avail¬ 
able  to  meet  future  requirements.  When  the  system  eventually  reaches 
steady  state,  we  find  that  we  now  have  excess  spares  on  hand  to  meet 
the  system  repair  requirements;  i.  e.  ,  five  failures  per  month. 

Although  the  new  power  supply  and  camera  units  were  similar 
to  the  old  system,  both  had  not  achieved  steady  state.  Therefore, 
assuming  that  our  previous  estimates  of  failure  rates  and  service 
times  were  applicable  to  the  new  systems  caused  us  to  acquire 
erroneous  stocks  of  parts/equipment  for  the  camera  unit. 

Totally  New  or  Dissimilar  Components 

When  the  Air  Force  procures  a  new  weapons  system,  a  large 
number  of  new  items  of  supply  are  also  received.  Since  the  Air  Force 
may  not  have  prior  supply  experience  for  these  new  units,  it  must 
depend  upon  the  contract  engineering  estimates  for  component  failure 
rates  and  service  times.  A  problem  arises  in  that  Air  Force  estima¬ 
tion  methods  may  not  be  the  same  as  those  used  by  the  contract 
engineers.  To  understand  this  problem,  we  need  to  define  a  few  key 
terms. 

If  we  are  discussing  a  population  of  components  that  operate 
together  as  a  system  in  steady  state,  and  if  each  component  is  repaired 
or  replaced  upon  failure,  we  use  the  term  Mean  Time  Detween 
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Failures  (MTBF).  (18:197)  "MTBF  may  be  determined  bytesting 
several  equipments  through  several  failure  incidents  and  averaging  the 
total  operating  time  for  the  number  of  failures.  "  (17:15-5)  If  we 
discuss  the  failure  or  service  time  of  one  component,  we  use  its  Mean 
Time  To  Failure  (MTTF).  We  could  life-test  several  of  these  units 
until  they  failed,  and  then  average  their  service  lives  (time  to  failure) 
to  find  the  unit's  MTTF. 

It  is  often  convenient  to  characterize  a  failure  model  or  set 
of  failure  data  by  a  single  parameter.  One  generally  uses 
the  mean  time  to  failure  or  the  mean  time  between  failures 
for  this  purpose.  .  .  .  Unfortunately  these  two  quantities  are 
sometimes  wrongly  thought  of  as  equivalent,  probably 
because  for  certain  simple  constant-hazard  cases  they  are 
equal.  In  a  single-parameter  distribution,  specification  of 
the  MTTF  fixes  the  parameter.  (18:197) 

In  Figures  1,  2,  and  3  the  power  supply  was  of  this  simple  type  where 

the  MTBF  calculated  by  our  steady  state  method  was  always  equal  to 

the  individual  unit's  MTTF.  This  was  not  true  of  the  camera  unit. 

When  procuring  a  new  system,  we  tell  the  contractor  exactly 

what  we  want.  We  want  a  component  that  will  have  a  given  MTBF 

(calculated  by  our  methods)  regardless  of  when  and  how  we  use  it. 

In  many  instances  military  contracts  for  complex  systems 
express  quantitative  reliability  requirements  in  terms  of 
mean  time  between  failures  (MTBF).  In  fact,  more  than  one 
MTBF  may  be  specified.  No  matter  how  many  MTBF's  are 
specified,  there  should  be  a  definition  of  what  MTBF  means. 

One  of  the  pitfalls  in  current  contracts  is  that  a  MTBF  is 
required  but  what  it  stands  for  is  not  explained.  .  .  . 

1.  Is  this  MTBF  the  design,  minimum  acceptable  hard¬ 
ware,  or  the  operational  requirement? 

2.  When  is  tt.is  MTBF  to  be  attained? 

3.  What  definition  of  failure  is  to  be  used  in  calculating 
MTBF? 
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The  second  question  above,  regarding  "when",  is  of  para¬ 
mount  interest.  The.  .  .  MTBF  may  be  attainable  at  some 
future  date,  but  not  in  the  first  operational  system. 

Usually,  the  best  approach  is  to  make  this  MTBF  the  aver¬ 
age  of  systems,  for  instance,  the  first  squadron  or  wing  of 
delivered  systems  for  a  military  contract.  .  .  .  (12:14-12) 

The  Air  Force  contracts  for  a  given  MTBF  and  the  contractor 
designs  for  a  steady  state  failure  rate.  But  how  long  will  it  take 
before  our  actual  failure  rate  reaches  the  steady  state  failure  rate? 

The  answer  to  this  may  be  different  for  each  component  population. 
However,  present  Air  Force  initial  provisioning  procedures  consider 
only  the  steady  state  MTBF.  Air  Force  planners  must  also  consider 
the  transition  period  .prior  to  steady  state.  To  do  this  we  must  know 
something  about  the  behavior  of  failure  rates  during  the  transition 
period. 

Burn-in  Phase 

There  is  another  problem  associated  with  the  present  "method.  " 
Many  complex  electronic  systems  go  through  an  infant  mortality  or 
"burn-in"  phase.  During  this  phase,  many  systems  exhibit  a  decreas¬ 
ing  failure  rate  with  respect  to  time.  Their  initial  failure  rate  is  high, 
but  after  the  first  few  repairs  the  system  failure  rate  settles  down  to  a 
constant  low  level.  (17:15-14)  This  condition  is  often  encountered  in 
the  real  world,  however  we  have  not  analyzed  this  type  of  phenomenon 
in  this  study.  A  long  run  average  or  an  engineer's  steady  state 
estimate  may  suppress  indications  of  a  burn-in  phase.  Unless  we  are 
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forewarned,  the  sudden  appearance  of  a  high  failure  rate  may  cause 
undue  alarm  during  the  first  few  operational  months. 

Summary 

In  this  chapter  we  have  attempted  to  show  that  the  procedures 
used  by  the  Air  Force  to  estimate  future  logistics  requirements  have 
necessary  underlying  assumptions.  The  most  basic  of  these  is  the 
assumption  of  a  steady  state  condition.  If  a  component  has  a  constant 
failure  rate,  its  component  poj  llation  will  immediately  appear  in 
steady  state.  However,  if  a  coi  lponent  has  an  increasing  or  a  decreas¬ 
ing  failure  rate  with  respect  to  time,  or  if  we  receive  a  steady  state 
engineering  estimate,  the  component's  population  may  pass  through  a 
transition  period  prior  to  reaching  the  steady  state  failure  rate.  The 
objective  of  this  research  is  to  investigate  failure  rates  to  be  expected 


during  this  transition  period. 
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General 

In  published  procedures  for  determining  quantitative  require¬ 
ments  for  initial  spare  parts,  the  Air  Force  assumes  that  the  occur¬ 
rence  of  component  failures  can  be  approximated  by  a  Poisson 
probability  distribution.  (5)  While  this  assumption  is  not  explicitly 
stated,  it  is  implicit  in  Air  Force  tables  which  suggest  that  every 
spare  part  has  a  unique  mean  time  between  failure  (MTBF)  and  a 
constant  failure  rate.  A  constant  failure  rate  implies  an  exponential 
density  function  (18:185)  and,  therefore,  that  failures  are  generated 
by  a  Poisson  process.  (8:111) 

In  the  following  discussion  we  will  show  that  there  are  two 
procedures  whi*_h  lead  to  a  Poisson  process.  The  first  procedure 
requires  each  component  to  have  a  constant  hazard  function,  but  does 
not  require  steady  state.  The  second  procedure  requires  a  steady 
state,  but  not  a  constant  hazard  function. 
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Constant  Hazard  Model 

If  we  are  discussing  one  component  and  we  can  assume  its  time 
to  failure  (TTF)  follows  an  exponential  distribution,  then, 

(1)  f(t)  =  Ae"^*  ,  0<  t  <  oo  (exponential  density  function) 

(2)  F(t)  =  f  f(x)dx  =  l-e"^*,  0  <  t  <  oo  (cumulative  density  function) 

*'0 


where  A  is  a  constant  and  is  defined  as  the  hazard  rate.  (18:48)  and  t 
is  operating  time.  V*e  can  also  find  the  mean  time  to  failure,  MTTF, 
by 


(3) 


MTTF  =  E 


[TTF]  =  f 


oo 

t  f(t)dt 


(7:165) 


However,  if  we  arc  discussing  a  population  of  identical  independent 
components,  we  can  discuss  the  fractional  portion  of  the  original 
population  which  fail  prior  to  a  time  t  by 


(4) 


P(x<  t)  =  f(x)dx  =  f 1  e_*Xdx  =  F (t), 
Jn  Jn 


<  t  < 


Which  is  identi  al  to  the  cumulative  function  above,  equation  (3).  Thus, 
1  -  F(t)  represents  the  fraction  of  the  population  surviving  to  time  t 
(7:164),  defined  as  the  reliability  function  as  follows, 

\S)  R(t)  =  1  -  F(t)  =  e"^1  (Reliability  function) 
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We  can  define  the  population  hazard  rate  or  the  distribution  hazard 
function  as  the  rate  at  which  the  survivors  fail.  (7 : 1 6 5 )  Then, 

(6)  Z{t)  =  =  A._— —  =  X  ,  0  <  t  <  oo  (hazard  function) 

R(t) 

which  is  constant  for  all  values  of  t.  For  this  distribution  it  has  been 
proven  that  the  reciprocal  of  the  hazard  function  is  the  mean  time 
between  failures,  (MTBF)  (14:4-43),  or, 

(?)  MTBF  =  — -  1 
Z(t)  X 

Thus,  for  the  constant  hazard  model  a  components  hazard  rate  equals 
the  population  hazard  rate  and  the  MTTF  =  MTBF,  and  the  MTBF  is 
independent  of  the  value  of  t. 

The  Renewal  Model 

If  we  are  discussing  the  failure  of  one  component,  but,  this 
time,  we  replace  the  failure  with  an  identical  new  component  (or  an  old 
component  which  has  been  rep?.iiv;u  to  new  condition),  we  have  a 
renewal  process.  (3)  Furthermore,  if  we  assume  time  to  be  component 
operating  time  and  that  we  can  replace  a  failed  component  instantan¬ 
eously,  the  system's  operating  time  becomes: 


<8)  Tsys  =  41  +  l2  +  4  3  *  "  +  fcn 


Where  tn  ~  the  operating  time  for  the  n^  renewal.  The  density  function 
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(p)  f(TSyS)  =  £(tl*t2-t3 . fcn) 

may  be  found  by  the  use  of  Laplace  transforms,  which  simplify  the 
function  to 

(10)  fTsys(S)  =  fl(S)  f2(S)  f3<S)  fn<S) 

where  f^ (S)  is  the  Laplace  transform  of  f n (t )  such  that 

(11)  fJ(S)  =^°°e-stn  fn(tn)dt 

or,  finally,  where  all  units  are  identical:^ 

(12)  4syg(S)  =  [f*(S)]n  (18:351-352) 

Up  to  this  point  we  have  been  discussing  the  renewal  (or  repair)  of  one 
component.  A  population  of  components,  all  operating,  fau  g,  and 
being  instantaneously  ’■enewed,  is  referred  to  as  a  superposition  of 
renewal  processes.  In  order  to  analyze  a  pool  of  individual  component 
processes,  it  is  helpful  to  visualize  each  component's  renewal  process 
as  an  independent  distribution  of  failure  times,  and  then  refer  to  the 
following  basic  reliability  references. 

^In  the  special  case  where  the  components  are  constant  hazard 
we  obtain: 

£n,  (t)  =  —  —  -  e  ^  (the  special  Erlangian  function) 

isys  (n-l) ! 


Which  is  called  the  Poisson  process. 
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One  of  the  first  areas  of  reliability  to  be  approached  with 
any  mathematical  sophistication  was  the  area  of  machine 
maintenance  (Khint chine,  1932,  and  C.  Palm,  1947).  The 
techniques  used  to  solve  these  problems  grew  out  of  the 
successful  experiences  of  A.  K.  Erlang,  C.  Palm  and  others 
in  solving  telephone  trunking  problems.  (2:1) 

As  Erlang  and  Palm  demonstrated,  the  Poisson  distribution  was  a 

very  good  approximation  to  the  actual  behavior  of  their  "simple 

stream"  of  phone  calls.  (11:12)  However,  Palm,  and  later,  Khintchine 

realized  that  groups  of  telephone  calls  did  not  always  satisfy  the  strict 

assumptions  of  the  Poisson  distribution.  This  was  especially  true  of 

the  requirement  that  the  distribution  be  memoryless  or  in  their  words, 

"without  after-effects.  "  (11:12)  In  the  introduction  to  his  mathematical 

proof  of  Palm's  theorem,  Khintchine  wrote: 

Thus,  although  an  investigator,  when  applying  the  theory 
to  experimental  data,  usually  has  the  task  of  explaining  the 
causes  of  deviation  of  the  actual  processes  from  the  course 
theoretically  predicted  for  them,  the  situation  is  in  fact  the 
wrong  way  round- -experimental  data  agree  with  the  deduc¬ 
tions  from  the  constructed  theory,  as  a  rule,  better  than 
could  have  been  expected  on  the  principal  considerations, 
and  it  is  precisely  this  'much  better'  agreement  which 
requires  explanation.  Palm  has  made  a  noteworthy  attempt 
to  explain  facts  of  this  kind  by  the  assumption  that  a  given 
process  represents  a  simple  sum  (superposition)  of  a  large 
number  of  mutually  independent  processes  o'  small  intensity, 
in  which  each  of  these  processes  is  stationary  and  orderly 
and  can  be  either  with  or  without  after-effects,  [author's 
underlining]  In  this  it  is  shown  that  under  broad  assump¬ 
tions  the  summary  stream  must  be  very  near  to  a  simple 
one.  Such  a  statement  of  the  problem  is,  on  the  face  of  it, 
very  near  to  the  real  situation.  So,  if  a  large  number  of 
subscribers  are  a  tached  to  the  given  exchange,  a  general 
stream  of  calls  is  composed  of  streams  (of  a  comparatively 
small  intensity)  arising  from  different  subscribers,  in  which 
the  component  streams  can  in  the  first  instance  be  regarded 
as  stationary,  orderly,  and  mutually  independent.  (11:49) 
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Khintchine  then  presents  a  chapter-long  proof  of  Palm's  proposal. 

While  this  particular  proof  is  too  lengthy  to  reproduce  in  this 

report,  two  other  authors  have  interpreted  the  meaning  for  us.  In  his 

text  on  renewal  theory,  Cox  discusses  Khintchine' s  proof  in  relation  to 

the  pooled  output  from  superpositioned  renewal  processes. 

Khintchine.  .  .has  proved  that;  in  the  limit  the  numbers  of 
renewals  in  non-cverlapping  intervals  follow  independent 
Poisson  distributions,  thus  showing  that  in  the  limit  the 
pooled  output  is  a  Poisson  process.  His  proof  does  not 
require  the  component  process  to  have  identical  distribu¬ 
tions  of  failure  time.  (3:78) 

And  in  his  text  on  reliability,  Shooman  writes: 

Since  we  are  now  dealing  with  the  superposition  of  sums 
of  random  variables,  the  arguments  for  using  the  central- 
limit  theorem  and  obtaining  limiting  forms  are  even 
stronger.  Specifically  one  can  state  that  for  .K  superposed 
(pooled)  processes: 

1.  If  each  process  is  Poisson  with  parameter  _X^,  the 
pooled  process  is  Poisson  with  parameter  K^  . 

2.  For  arbitrary  processes,  if  a  large  number  is  pooled, 
the  output  approaches  a  Poisson  process. 

3.  As  K  becomes  large,  the  time  between  renewals 
approaches  an  exponential  distribution  with  parameter 
K\.  .  .  (18:356-357) 

If  each  component  of  a  population  operates,  fails,  and  is  renewed;  we 
have  a  population  of  independent  renewal  processes.  If  such  a  popula¬ 
tion  were  in  the  Air  Force  inventory,  the  Air  Force  would  pool  the 
outputs  from  these  processes  into  a  maintenance  report.  Khintv  .ne's 
and  Shooman' s  proofs  show  that  this  procedure  would,  in  the  limit, 


produce  a  Poisson  process. 
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Summary 

We  have  attempted  to  demonstrate  in  this  chapter  that  there  are 
two  methods  which  produce  a  Poisson  process.  First,  if  a  component 
has  a  constant  failure  rate,  it  can  be  described  by  the  constant  hazard 
model.  This  model  is  independent  of  time  and  number  in  the  popula¬ 
tion.  Second,  if  a  component  has  an  increasing  or  decreasing  failure 
rate  with  time,  the  constant  hazard  model  does  not  apply.  However, 
the  Poisson  process  still  applies,  at  steady  state,  if  the  component  is 
part  of  a  pooled  renewal  process. 


CHAPTER  IV 


MATHEMATICAL  TOOLS 


The  Weibull  Distribution 


The  Poisson  and  exponential  distributions  enjoyed  wide  accept¬ 
ance  in  the  1950's  due  to  their  good  results  and  appealing  simplicity. 

A  comprehensive  analysis  of  failure  distributions  was  done 
[in  the  early  1 950 '  s]  at  [the]  RAND  [Corporation]  by  D.  J. 

Davis,  who  found  that  the  exponential  distribution.  .  .accur¬ 
ately  described  the  failure  characteristics  of  a  v/ide  variety 
of  devices.  As  examples,  he  cited.  .  .  commercial  aircraft 
radio  tubes,  radar  set  components,  tubes,  resistors,  .  .  . 
and  combinations  of  various  components.  (9:209) 

Many  other  investigators  added  items  to  the  growing  list.  Yet,  there 

were  intuitive  and  theoretical  arguments  against  the  universal  use  cf 

these  distributions. 

The  exponential  law  is  reasonably  appropriate  where  chance 
alone  dictates  failure  occurrence.  If  it  is  known,  for  example, 
that  failure  is  consistently  due  to  deterioration,  wear-out, 
degradation,  fatigue,  or  any  Repetitive  mechanism,  almost 
certain  m.noxoojientiality  is  implied.  .  .  .  To  be  truly  random, 
failure  cannot  be  due  to  design  deficiencies  or  manufacturing 
errors.  Being  assignable  as  to  cause,  failure  distributions 
associated  with  early  life  of  a  product  usually  are  not  expo¬ 
nential  in  that  they  are  not  random.  (13:3^4) 

O. 

We  were  cautioned -- 
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It  is  to  be  emphasized  that  considerable  care  is  required 
in  the  selection  of  the  appropriate  underlying  distribution 
for  reliability  testing.  The  validity  of  the  test  results 
depends  to  a  large  degree  upon  how  we.’l  the  selected  prob¬ 
ability  distribution  represents  the  actual  distribution  of  the 
time  to  failure  upon  which  observations  are  being  made.  .  .  . 

This  point  is  emphasized  because  the  widespread  (perhaps 
indiscriminate)  use  of  the  exponential  distribution  as  a  model 
of  failure  patterns  may  lead  one  to  believe  that  failure  times 
in  general  may  be  adequately  represented  by  such  a  distri¬ 
bution.  .  .  .  The  use  of  other  distributions  introduced  compli¬ 
cations  into  the  testing  procedures;  however,  if  one  is  to 
draw  reasonable  conclusions  from  the  tests,  he  must  he 
willing  to  face  these  other  problems.  (4:3-6) 

Professor  Weibull  was  willing  to  face  the  problems.  He  ".  .  . 
suggested  it  [the  Weibull  distribution]  as  the  distribution  for  a  variety 
of  problems.  .  .  .  The  agreement  he  found  between  his  observations  and 
those  predicted  with  the  fitted  Weibull  was  very  impressive.  ..."  (1:1) 

The  Weibull  Mathematical  functions  are: 

Weibull  density  function: 


(13)  f(f;)  =  JL  (t-c)k-1  EXP 

0k 


[-It]'] 


,  K>O,0>O,O  <  t< 


00. 


Weibull  cumulative  function: 


(14)  F(t)  =  f  f(t)dt  =  1  -  EXP 

J  r\ 


■  [t-d 

k 

.  u 

,  K>O,0>O,O<t<  oo. 


Weibull  reliability  function: 


(15)  R(t)  =  EXP 


r  t~c  k 

Q 


,  K  >  0,  0>  o ,  0  <  t  <  00. 


Where  C  is  the  location  parameter  (for  purposes  of  this  research 
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C  =  0  and  time  begins  at  the  origin);  K  is  the  shape  parameter;  and, 

6  is  the  scale  parameter. 

Since  this  distribution  includes  the. .  .  exponential  as  a  special 
case,  the  Weibull  is  sometimes  thought  of  as  a  generalization 
of  the  exponential  distribution.  Its  use  would  then  be  consid¬ 
ered  by  anyone  who  has  been  using  the  exponential  distribu¬ 
tion  but  felt  a  more  flexible  model  was  needed.  Indeed,  the 
popularity  of  the  exponential  distribution  as  the  model  for 
many  experimenters  in  the  fifties  was  challenged  by  the 
popularity  of  the  Weibull  in  the  sixties.  (1 : Iff) 

The  Weibull  distribution,  with  its  generality,  has  indeed  given 
the  model  builder  a  powerful  tool.  By  the  appropriate  choice  of  the 
parameter,  K,  we  can  change  the  distribution  from  the  decreasing 
hazard  model--a  model  sometimes  used  to  show  the  "infant  mortality" 
or  break-in  period,  to  the  constant  hazard  exponential  model,  to  the 
increasing  hazard  model--models  characterized  by  the  Rayleigh  and 
normal  distribution  where  a  unit  wears-out  with  time.  By  proper 
choice  of  the  scaling  factor  6,  we  can  adjust  the  amplitude  of  the  dis¬ 
tribution  much  like  we  do  when  choosing  the  appropriate  X  for  the 
exponential.  Therefore,  by  adjusting  K  and  6  we  may  model  a  great 
number  of  hazard  distributions.  (18:190) 

To  properly  define  the  failure  rate  characteristics  of  a  system 
of  equipment/ components,  it  is  necessary  to  describe  the  instantaneous 
failure  rate  function  associated  with  the  underlying  failure  generating 
distribution.  This  instantaneous  failure  rate  function  has  been  gen¬ 
erally  described  by  many  authors  as  the  hazard  rate  function,  Z(x), 


as  follows: 


(16) 


(10:2-5) 
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Where  f(x)  is  the  failure  density  function  and  R(x)  is  the  reliability 
function. 

The  failure  density  function  has  been  defined  as: 

(17)  f(x)  =  Z(x)  EXP  Z(x)dxj  (10:2-5) 

Thus  given  an  instantaneous  failure  rate  function,  we  can  easily  obtain 
the  failure  density  function  by  simple  integration  based  upon  the 
exponential  law  of  growth  or  decay. 

Many  failure  distributions  can  therefore  be  determined  by 
integration  of  their  characteristic  instantaneous  failure  rate  function. 
For  example: 

If:  Z(t)  =  |*0>O 

Then:  f(t)  =  I  EXP  [-  f  I  dx 

*  lJo  e  J 


which  is  the  well  known  exponential  function. 

Other  instantaneous  failure  rate  functions  may  be  integrated  to 
obtain  the  required  failure  density  functions  in  the  same  manner. 

The  instantaneous  failure  rate  function  used  in  this  simulation 


is  as  follows: 
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(18)  Z(t)  '=  -T-  tk-1,  $>  0,  K>  0, 0  <  t  <  ”, 

eK 


Then,  f(t)  =  Z(t)  EXP 


Xt  1 

Z(t)dt 

u 


(19)  f(t)  =  H-  tk_1  EXP 

0k 


j  »  0  >  0  *  ^>0 


,  0  <  t  <  03  . 


This  is  one  form  of  the  Weibull  function  and  it  is  interesting  to  note 
that  the  exponential  distribution  is  a  special  case  of  the  Weibull  distri¬ 
bution  when  K  =  1.  Other  authors  have  indicated  that  when  K  is 
approximately  3,  the  Weibull  approximates  the  shape  of  the  Normal 
distribution.  (19:21)  Other  failure  distribution  shapes  can  be  approxi¬ 
mated  by  the  Weibull  model.  Decreasing  failure  rate  models  have  a 

0  <  K  <  1  and  increasing  failure  rate  models  have  a  K  >  1. 

The  parameter  6  of  the  Weibull  model  determines  the  variance 
or  dispersion  of  the  model.  The  variance  of  the  Weibull  model  is 
given  as: 

(20)  Variance  Weibull  =  02jr(l  +  2/k)  -  (T(l  +  1/k))2 

Thus  by  changing  the  Rvalue  we  can  set  the  variance  of  the  failure 
density  function  and  in  conjunction  with  the  K  value  approximate  many 
failure  distributions. 


Mathematics  of  the  Model 

Using  equations  16,  17,  and  18  of  the  previous  section,  we  can 
derive  the  following  mathematical  tools  used  in  this  study. 
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Hazard  function: 


(21)  Z(t)  =  |k  tk_1,  K  >0  ,6>  0,  0  <  t  < 


00* 


Failure  density  function: 


(22)  f(t)  =  Z(t)  EXP 


^tk  1  EXP 


ft)k' 

PI 


,  K>O,0>O,O<t<oo. 


Cumulative  density  function: 


/»t 

(23)  F(t)  =  /  f(t)dt,  0<  t<  «. 
Jc\ 


=  1-EXP 


•MY 


,  K>0,  0>O,  0  <  t  <  oo 


Reliability  function: 


(24)  R(t)  =  1  -  F(t) 


EXP  jj- (t / (9  )kj  ,  K>0,  6>0,  0  <  t  <  oo. 


where  K  is  the  shape  parameter  and  6  is  the  scale  parameter.  (15:94) 
From  equation  3  and  22  we  can  then  determine  the  expected  value  of 
the  density  function,  MTTF,  as  follows: 


Jr  00 

t  f(t)dt 

o 


r  00 ,  K  .k-i  r 

■jC*  Exp 

« 


Idt 


=  er(  1  +  1/k),  0>O,  K  >  0 


As  shown  earlier,  this  relationship  yields  the  steady  state  failure  rate 
for  a  pooled  output  of  renewal  processes,  where  each  component's 
density  function  is  described  by  a  Weibull  density  function  with 
parameters  K  and  6 ,  and  if  the  number  of  processes  is  large.  Under 
such  conditions,  we  can  equate  the  MTTF  of  the  individual  components 
to  the  MTBF  of  the  component  population  as  follows: 

(25)  MTBF  =  dr  (1  +  1/k) 


We  used  equation  25  to  compute  our  simulation  model  scale 
parameter,  6.  When  the  shape  parameters,  0.  5  5  K<  5.  0,  were  then 
chosen,  we  demonstrated  that  an  appropriate  choice  of  the  scale 
parameter,  Q,  caused  the  pooled  output  of  the  renewal  processes  to 
approach  the  MTBF  predicted  by  the  constant  hazard  function  of  the 
exponential,  :,e., 

MTBFgxp  -  Constant  Failure  Rate  * 

And  therefore,  at  steady  state,  the  MTTF  for  any  individual  component 
would  equal  the  MTBF  for  the  component  population--a  condition  which 
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exists  only  in  a  Poisson  process.  (18:197) 


Random  Numbers 

To  generate  random  numbers  to  compute  individual  times  to 
failure  for  components,  we  used  the  inverse  transform  method  as 
follows:  (20:511) 


v  =  1  -  EXP 


EXP 


1-v  =  u 


-  (t/0)k  =  In  u 


tk  =  0k(-ln  u) 


t  =  6  (-In  u)1/k 


where  v  is  a  uniformly  distributed  pseudo-random  number,  0  <  v  <  1, 
and  thus,  (1  -  v)  =  u  is  also  a  uniformly  distributed  pseudo-random 
number,  0  <  u  <  l. 

Therefore,  by  this  inverse  transform  procedure,  we  can 
generate  a  stream  of  values  for  t  which  follow  the  Weibull  distribution. 


mmmm  m  m 


mmum* 


l|.i  iKWIPPViil  M 


CHAPTER  V 


PROCEDURES 


The  Computer  Model 

We  developed  a  computer  model  for  this  investigation  to  simu¬ 
late  the  renewal  processes  underlying  a  portion  of  the  Air  Force 
maintenance  system. 

The  time  frame  of  our  model  covers  approximately  four  years. 
We  selected  a  30  day  month  as  our  reporting  interval  so  that  our  4 
year  period  contained  48  monthly  reporting  points.  We  assumed  v.  4 
hour  daily  utilization  rate  per  component  and,  thus,  each  component  in 
the  inventory  operated  120  hours  each  month.  The  required  inputs  for 
each  simulation  were  the  total  fleet  size,  the  introduction  rate,  and 
the  desired  shape  parameter. 

The  computer  model  introduced  a  percentage  of  the  total  fleet 
into  our  inventory  at  the  beginning  of  each  month.  As  each  unit 
entered,  a  service  time  was  selected  from  a  random  Weibull  distribu¬ 
tion  and  then  added  to  the  present  computer  time  to  form  the  unit's 
future  failure  time.  The  nodcl  then  stored  the  unit  with  its  failure 
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time  in  a  future  events  file.  When  the  computer  clock  reached  a 
scheduled  failure  time,  the  model  removed  the  associated  unit  from 
the  events  file;  recorded  the  failure;  selected  a  new  future  failure  time; 
and  refiled  the  unit  and  time  in  the  events  file.  Once  started,  this 
cycle  continued  independently  for  each  unit  in  the  inventory. 

At  the  end  cf  each  computer  month,  the  model  stored  the  total 
number  of  failures  which  occurred  in  a  memory  array.  We  repeated 
each  4  year  simulation  40  times  for  each  set  of  input  parameters. 

After  the  40th  run,  the  model  calculated  an  average  number  of  failures 
for  each  monthly  period,  a  id  then  divided  this  value  by  each  month's 
operating  time.  The  results  were  48  average  monthly  failure  rates. 

For  a  more  detailed  description  of  the  simulation  model  see  Appendix  A. 

Scope 

General.  --The  research  for  this  paper  was  limited  to  a  litera¬ 
ture  search  and  digital  computer  simulation  techniques.  We  have  not 
attempted  to  evaluate  our  simulated  data  against  actual  failure  data. 
However,  we  have  attempted  to  select  our  model  parameters  so  that 
they  would  present  a  reasonable  approximation  to  actual  conditions. 

Fleet  Size.  -  -The  total  size  of  a  fleet  of  aircraft  may  depend  on 
many  variables,  e.  g. ,  planned  mission,  projected  costs,  actual  cost, 
adaptability  to  alternate  missions,  world  conditions,  etc.  Therefore, 
we  set  a  range  of  tleet  sizes,  20  to  1000  aircraft,  and  tested  incre¬ 
ments  of  this  range  in  our  preliminary  simulations.  Our  preliminary 
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results  indicated  that  the  total  number  in  the  fleet  had  little  effect  on 
our  prime  performance  variable,  failure  rate,  as  long  as  a  constant 
percentage  of  that  total  were  introduced  each  month.  For  this  reason 
we  chose  a  fleet  size,  100  aircraft,  that  was  well  within  the  feasible 
range,  and  one  that  would  simplify  later  extrapolation. 

Introduction  Rate.  -  -As  with  fleet  size,  the  introduction  rate  of 
a  fleet  of  aircraft  is  subject  to  political  and  economic  pressures  and, 
therefore,  highly  variable.  In  addition,  this  rate  is  limited  by  the 
physical  problems  of  "tooling  up"  for  production.  Therefore,  we 
selected  as  an  average  figure  the  monthly  introduction  rate  at  5%.  We 
used  percentage  values,  as  stated  before,  to  simplify  extrapolation  to 
actual  introduction  phases. 

Average  Operating  Period.  -  -The  number  of  hours  ?.n  aircraft 
flies  each  day  again  depends  on  the  type  of  aircraft  and  the  mission. 
However,  when  we  are  discussing  components  on  an  adrcraft,  the 
problem  becomes  more  complex.  A  C-5A  may  use  a  system  for  the 
full  mission,  a  B-52  may  use  the  same  type  system  only  on  take  off 
and  landing,  and  an  F-106  may  use  the  same  system  for  special 
missions  only.  In  our  model,  all  component  failures  are  a  function  of 
operating  time;  thus,  all  results  are  directly  related  to  the  utilization 
rate.  For  this  reason  we  chose  the  "ball  park"  value  of  four  hours 
per  unit  per  day  for  our  utilization  rate. 

Mean  Time  Between  Failures.  -  -We  selected  the  values  for  this 
variable  to  cover  a  feasible  range  for  Air  Force  repair  or  replacement 
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parts.  We  conjectured  that  a  major  system  with  less  than  50  hours 
MTBF  would  be  unacceptable  for  reliability  reasons,  and  that  a  system 
with  failure  intervals  greater  than  2000  hours  (nearly  17  months  in 
our  model)  would  cause  few  logistics  problems.  We  then  investigated 
the  values  of  50,  100,  500,  1000,  and  2000  hours  to  cover  our  assumed 
range. 

Shape  Parameters.  --We  investigated  the  effects  of  this 
parameter  through  the  values  of  0.  5,  1.0,  1.5,  2.0,  3.0,  4.0,  and 
5.  0;  and  we  conjecture  that  the  properties  of  our  performance  variable, 
failure  rate,  are  continuous  within  this  range.  These  values  take  us 
from  a  rapidly  decreasing  failure  rate  (infant  mortality  phase)  through 
a  constant  failure  rate  (random  failures)  to  a  rapidly  increasing  failure 
rate  (wear-out  phase). 

Scale  Parameters.  --The  scale  parameter,  $,  was  determined 
as  previously  indicated  from  equation  25.  Equation  25  allowed  us  to 
transform  the  MTBF  desired  or  indicated  to  an  appropriate  scale 
parameter  of  our  simulation  model.  Thus,  we  were  able  to  track  the 
experienced  failure  rates  through  the  transition  period  until  the  com¬ 
ponent  population  steady  state  failure  rate  was  achieved. 


CHAPTER  VI 


SIMULATION  RESULTS 


General 

The  simulation  results  have  been  divided  into  decreasing 
hazard  rate,  constant  hazard  rate  and  increasing  hazard  rate  cate¬ 
gories  to  more  easily  discuss  and  analyze  the  simulation  output.  The 
discussion  and  analysis  will  be  based  upon  frequency  plots  of  the  out¬ 
put  data.  The  actual  tabulated  output  data  is  included  in  Appendix  B 
for  readers  wishing  to  perform  a  more  extensive  analysis  than  p  re¬ 
sented  in  this  study. 

As  indicated  earlier,  our  objective  is  to  investigate  the  distri¬ 
bution  of  failures  per  hour  during  the  transition  period  from  initial 
acquisition  of  weapons  systems  to  the  attainment  of  steady  state;  and 
we  want  to  compare  such  failure  rates  to  the  constant  failure  rate 
predictions  presently  in  use  by  the  Air  Force.  (5)  Thus,  our  main 
interest  is  the  magnitude  of  the  difference  between  the  simulation 
results  and  the  Air  Force  prediction  methods. 
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Definitions 

The  analysis  has  been  performed  in  terms  of  the  average  abso¬ 
lute  deviation  between  the  present  prediction  methods  anc  the  simula¬ 
tion  results  prior  to  the  achievement  of  the  system  steady  state  condi¬ 
tion.  For  purposes  of  analysis,  system  steady  state  is  "broadly" 
defined  as  relatively  minor  absolute  deviations  of  the  simulation 
results  from  the  current  predicted  constant  failure  rate. 

The  following  definitions  are  necessary  to  the  analysis. 

Constant  Hazard  Rate  (CHR)--the  constant  failure  rate  predicted 
by  present  Air  Force  methods 

CHR  = _ 1 _  . 

MTB^ 


Constant  Failure  Rate- -the  condition  of  the  Weibull  failure 
generating  function  where  the  shape  parameter,  K,  is  equal  to  1. 

Decreasing  Failure  Rate--the  condition  ot  the  Weibull  failure 
generating  function  where  the  shape  parameter,  K,  is  0  <  K<  1. 

Increasing  Failure  Rate--the  condition  of  the  Weibull  failure 
generating  function  where  the  shape  parameter,  K,  is  greater  than  1. 

Mean  Time  Between  Failure  (MTBF)--the  time  between  failures 


predicted  by  the  present  Air  Force  method 


MTBF  =  — 


Constant  Failure  Rate 
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average  absolute  deviation  between  the  simulation  results  and  the  CHR 
computed  as  follows; 

48 

T  |  PFR:  -  CHR  I 

AD  =  ifi _  x  100%  . 

48  (CHR) 

Per  Cent  Maximum  Absolute  Deviation  (MD)--the  per  cent  of 
maximum  absolute  deviation  between  the  simulation  results  and  the 
CHR  computed  as  follows: 

MD  =  MAX]  PFRi  -  CHR]  x  100%,  i  =  1,  2,  3,  .  .  .  ,  48. 

CHR 

Per  Cent  In  Range  (IR)--the  number  of  indivi  »aal  PFR's  within 
the  range  CHR  +  AD  computed  as  follows: 

48 

Z  Oi 

IR  =  ill -  x  100o,o 

48 

where: 

[  1,  when  (CHR  -  AD  )<PFR.  <  (CHR  +  AD) 

0,  otherwise 

The  Decreasing  Failure  Rate  Case 

The  Decreasing  Failure  Rate  Case  is  defined  as  the  condition 
where  the  hazard  function  of  the  components  is  inversely  related  to 
time,  i.  e.  ,  the  function  decreases  over  time.  In  this  case,  the  simu¬ 
lation  program  assigned  to  each  component  of  the  population  such  a 
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decreasing  failure  rate,  and  we  were  interested  in  observing  what 
effect  such  a  failure  rate  would  have  on  the  total  population  failure 
rate. 

Figures  6,  7,  and  8  of  Appendix  B  are  plots  of  the  system 
failure  rates  produced  by  the  simulation  program.  (The  simulation 
output  may  be  found  in  Appendix  B  of  this  study.  )  The  plots  of  the 
system  failure  rates  indicate  relatively  large  failure  rates  during  the 
initial  months  of  introduction.  However,  these  large  failure  rates 
rapidly  reduced  in  magnitude  in  a  relatively  short  period  of  time.  The 
rate  of  change  of  the  system  failure  rate  is  negative  and  decreases  at 
a  decreasing  rate  eventually  approaching  the  CUR  asymptotically  from 
above. 

Each  plot  indicates  that  deviations  from  the  CHR  may  be  10% 
or  higher  for  the  first  24-28  months  of  operation.  Such  deviations 
appear  to  be  directly  related  to  the  size  of  the  MTBF  as  can  be  seen 
from  the  plots.  Figure  6  indicates  that  with  a  MTBF  of  500  hours,  the 
system  of  components  enters  steady  state  (our  broad  definition)  in 
approximately  36  months.  However  Figures7  and  8  with  MTBF's  of 
1000  hours  and  2000  hours  are  not  to  be  considered  in  steady  state 
throughout  the  entire  simulation  period  of  48  months. 

The  simulation  results  indicate  that  components  of  a  given 
population  characterized  by  a  decreasing  hazard  function  affect  the 
population  failure  rate  distribution.  Much  higher  population  rates  than 
predicted  by  the  present  Air  Force  methods  should  be  experienced 
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during  the  initial  months  of  the  transition  period  to  steady  state.  It 
appears  that  the  population  failure  rate  will  eventually  achieve  steady 
state  conditions,  however,  the  transition  period  may  be  quite  lengthy, 
and  for  large  MTBF's  may  not  actually  be  attainable.  Table  1  indicates 
the  simulated  maximum  %  deviations  (MD)  and  the  average  %  deviations 
(AD)  from  such  a  CHR  prediction  method. 


TABLE  1 

THE  DECREASING  FAILURE  RATE  CASE 


MTBF 

Shape 

Parameter 

MD 

AD 

%  OBS. 

(1/MTBF)  (1+AD) 

500 

.5 

272.  9% 

36.57% 

68.  75% 

.5 

412.  5% 

59.89% 

60.41% 

2000 

.5 

591.6% 

95.58% 

62.50% 

Under  such  conditions,  the  Air  Force  planners  may  decide  to 
approximate  the  expected  failure  rate  by  averaging  several  period 
failure  rates.  However,  since  the  population  failure  rate  distribution 
is  continuously  decreasing,  such  a  procedure  would  result  in  the 
acquisition  of  replaccment/spare  parts  in  quantities  larger  than 
actually  necessary.  If  the  MTBF  supplied  by  the  contractor  is  used 
as  the  basis  for  logistical  support  planning,  the  simulation  results 
indicate  that  fewer  items  would  be  planned  for  than  would  actually  be 


required. 
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The  Constant  Failure  Rate  Case 

The  Constant  Failure  Rate  Case  is  defined  as  the  condition 
where  the  hazard  rate  function  of  the  components  is  a  constant  over 
time.  The  simulation  program  assigned  to  each  component  of  the 
population  such  a  failure  rate.  Again,  the  effect  of  such  a  component 
characteristic  on  the  total  population  failure  rate  is  our  main  interest. 

Figures  9,  10,  11,  12,  and  13  of  Appendix  B  are  plots  of  the 
system  failure  rates  produced  by  the  simulation  program.  Al.’  plots 
indicate  that  the  system  failure  rates  immediately  tend  to  oscillate 
around  the  CHR  predicted  by  the  present  Air  Force  method.  A  transi¬ 
tion  period  does  not  appear  to  be  present  and  according  to  our  "broad" 
definition  of  steady  state,  we  conjecture  the  system  of  components  to 
be  immediately  in  steady  state.  The  oscillations  around  the  predicted 
CHR  tend  to  increase  proportionally  with  increases  in  the  MTBF  used. 
Table  2  indicates  the  simulated  maximum  %  deviations  (MD)  and  the 
average  %  deviations  (AD)  from  the  CHR. 


TABLE  2 

THE  CONSTANT  FAILURE  RATE  CASE 


MTBF 

Shape 

Parameter 

MD 

AD 

%  OBS 

(1/MTBF)  (1+AD) 

50 

1.  0 

8.  33% 

1.  78% 

56.  25% 

100 

1.  0 

10.  83% 

2.  06% 

62.  50% 

500 

1.0 

10.95% 

3.  14% 

64.  58% 

1000 

1.  0 

17.40% 

4.  79% 

58.  33% 

2000 

1.0 

20.  80% 

5.48% 

60.41% 
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Under  such  conditions  it  appears  that  the  present  Air  Force 
procedures  for  predicting  expt  ted  failures  for  new  weapons  systems 
would  produce  "satisfactory"  results  since  the  system  would  immed¬ 
iately  be  in  a  steady  state  condition. 

The  Increasing  Failure  Rate  Case 

The  Increasing  Failure  Rate  Case  is  defined  as  the  condition 
where  the  hazard  function  of  the  components  is  directly  proportional 
to  time,  i.  e. ,  the  function  increases  over  time.  As  with  the  other 
cases,  we  are  still  interested  in  the  effect  such  component  failure 
rates  have  upon  the  population  failure  rate.  We  conjecture  that  com¬ 
ponent  failures  which  possess  this  characteristic  are  more  represent¬ 
ative  of  the  real  world  situation  than  either  the  decreasing  or  constant 
failure  rate  cases.  Intuitively,  it  i‘  more  reasonable  to  expect  larger 
probabilities  of  failure  as  items  age  tt:?n  either  constant  or  decreasing 
probabilities  of  failure  with  increasii  g  age. 

Figures  14  through  36  of  Appendix  B  are  plcts  of  the  simulation 
results  fcr  varying  MTBF's  and  shape  parameters  for  this  case.  (The 
simulation  data  may  be  found  in  Appendix  B  of  this  study.  )  The  plots 
of  the  simulated  failure  rates  have  been  categorized  into  groups  on  the 
basis  of  MTBF’s  and  will  be  discussed  in  this  order. 

MTBF  =  50.  --For  shape  parameters,  2<k<5,  the  plots  indi¬ 
cate  a  minimum  of  5%  deviation  between  the  predicted  failure  rate  and 
the  simulation  results  for  the  first  few  months  during  transition  to 
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steady  state.  The  magnitude  of  these  deviations  appears  to  increase 
with  increases  in  the  shape  parameter.  Such  a  condition  is  to  be 
expected  with  the  Weibull  distribution,  since  shape  parameters  less 
than,  or  greater  than  approximately  3  generate  skewed  distributions, 
right  and  left  respectively,  with  a  symetrical  distribution  occurring 
when  the  shape  parameter  equals  3.  The  length  of  the  transition  period 
appears  to  be  in  the  16-24  month  range  for  MTBF's  of  50,  and  the 
variation  around  the  predicted  CHR  during  steady  state  is  small.  The 
rate  of  change  in  the  population  failure  rate  distribution  is  positive  and 
approaches  the  predicted  CHR  asymptotically  from  below.  See  Figures 
14  through  17,  Appendix  B. 

MTBF  =  100. --Again  the  same  general  observations  can  be 
made  of  these  plots.  The  rate  of  change  of  the  failure  rate  is  positive, 
increasing  at  a  decreasing  rate  and  approaching  the  predicted  CHR 
asymptotically  from  below.  The  system  failure  rate  appears  to  achieve 
steady  state  in  about  20  months.  The  simulated  results  deviate  from 
the  predicted  CHR  by  at  least  10%  during  the  first  few  months  of  the 
transition  period.  The  maximum  deviation  of  the  simulated  failure  rate 
from  the  predicted  CHR  is  approximately  31%  during  the  first  month. 
See  Figures  18  through  21,  Appendix  B. 

MTBF  =  500.  --The  same  general  observations  made  above  are 
again  observed.  However,  the  maximum  deviation  of  the  simulated 
failure  rates  from  the  predicted  CHR  has  increased  considerably.  For 
shape  parameters  of  4  and  5  the  maximum  deviation  is  100%.  The 
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system  appears  to  achieve  steady  state  in  approximately  24-28  months. 
The  rate  of  change  in  the  failure  rate  also  appears  to  have  decreased 
with  increasing  scale  parameters.  See  Figures  22  through  26, 

Appendix  B. 

MTBP  =  1000.  The  same  general  observations  made  above  also 
apply  here.  Steady  state  appears  to  be  achieved  in  approximately  24-28 
months,  however,  the  relative  deviations  around  the  predicted  CHR  at 
steady  state  has  increased.  The  rate  of  change  of  the  failure  rate  has 
again  decreased.  See  Figures  27  through  31,  Appendix  B. 

MTBF  =  2000.  The  same  general  observations  are  applicable 
to  this  group  also.  There  is  a  marked  difference  in  the  rate  of  change 
of  the  system  failure  rate.  Steady  state  appears  to  be  reached  in  the 
26-28  month  period.  The  relative  fluctuations  in  the  failure  rate  have 
also  increased  markedly.  During  the  first  12  months,  the  minimum 
deviation  between  the  simulation  results  and  the  predicted  CHR  is  at 
least  40%.  See  Figures  32  through  36,  Appendix  B. 

In  all  of  the  plots,  there  appears  to  be  a  change  to  the  shape  of 
the  failure  distribution  during  the  transition  period  resulting  from 
increasing  MTBF's.  The  relative  fluctuation  of  simulated  failure  rates 
also  appears  to  increase  as  the  MTBF  is  increased.  All  the  failure 
distributions  approach  t'  predicted  CHR  asymptotically  from  below. 
Tables  3  and  4  summa  ze  the  results  of  the  simulation  for  the 
parameters  and  ranges  tested.  Table  3  indicates  that  for  increases  in 
MTBF  for  a  given  shape  parameter,  the  MD  and  AD  (as  defined  before; 


46 


TABLE  3 

INCREASING  FAILURE  RATE  CONDITION  WITH  CONSTANT  SHAPE 


MTBF 

Shape 

Parameter 

MD 

AD 

%  OBS 

(1/MTBF)  (1+AD) 

50 

2.  0 

12.  71% 

1.43% 

77.  08% 

100 

2.  0 

25.  42% 

2.  73% 

75.00% 

500 

2.  0 

77.  10% 

10.  63% 

70.  80% 

1000 

2.  0 

85.  40% 

18.  35% 

60.41% 

2000 

2.  0 

100.00% 

28.  74% 

60.41% 

50 

3.0 

15.83% 

1.  54% 

79.  16% 

100 

3.  0 

31 . 67% 

2.  94% 

79.16% 

500 

3.  0 

93.  75% 

13.96% 

72.91% 

1000 

3.  0 

100.00% 

21. 7 1 % 

62.  50% 

2000 

3.  0 

100.00% 

34.  46% 

62.  50% 

50 

4.  0 

16.  66% 

1.  56% 

75.00% 

100 

4.0 

31.67% 

2.  94% 

72.91% 

500 

4.0 

100.00% 

14. 47% 

70.83% 

10"0 

4.  0 

100.00% 

23.  09% 

64.  58% 

2000 

4.  0 

100.00% 

36.44% 

62.  50% 

50 

5.  0 

17.5% 

1.  58% 

77.  08% 

100 

5.  0 

30.00% 

2.99% 

72.91% 

500 

5.0 

100. 00% 

14.  51% 

72,91% 

1000 

5.  0 

100.00% 

24. 06% 

64.  58% 

2000 

5.0 

100.00% 

37.  70% 

62.  50% 

have  increased.  Table  4  indicates  that  for  increases  in  the  shape 
parameter  for  a  given  MTBF  that  the  MD  and  AD  have  increased.  Thus 
for  large  MTBF's  or  large  shape  parameters,  we  can  expect  larger 


deviations  from  the  predicted  CHR. 
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TABLE  4 

INCREASING  FAILURE  RATE  CASE  WITH  CONSTANT  MTBF 


MTBF 

Shape 

Parameter 

MD 

AD 

%  OBS 

(1  /MTBF)  (1+AD) 

50 

2.0 

12.71% 

1.43% 

77.  08% 

50 

3.  0 

15.83% 

1.  54% 

79. 16% 

50 

4.  0 

16.  66% 

1.  56% 

75.  00% 

50 

5.  0 

17.  50% 

1.58% 

77.  08% 

100 

2.  0 

25.42% 

2.  74% 

75.  00% 

100 

3.  0 

31.  67% 

2.95% 

79.  16% 

100 

4.  0 

31.  67% 

2.  94% 

72.91% 

100 

5.  0 

30.  00% 

2.  99% 

72.  91% 

500 

2.  0 

77.  10% 

10.  63% 

70. 80% 

500 

3.  0 

93.  75% 

13.96% 

72.91% 

500 

4.  0 

100. 00% 

14.  47% 

70.  83% 

500 

5.  0 

100.00% 

14.51% 

72.91% 

1000 

2.  0 

85.40% 

18.  35% 

60.41% 

1000 

3.  0 

100.00% 

21.71% 

62.50% 

1000 

4.  0 

100.00% 

23.  09% 

64.  58% 

1000 

5.  0 

100.  00% 

24.  06% 

64. 58% 

2000 

2.  0 

100.00% 

28.  74% 

60.41% 

2000 

3.  0 

100.00% 

34.  46% 

62.  50% 

2000 

4.  0 

100. 00% 

36.44% 

62.  50% 

2000 

5.  0 

100. 00% 

37.  70% 

62.  50% 

Standardised  Table  for  Steady  State 

The  following  is  the  development  of  a  standard  table,  Table  5, 
and  curve,  Graph  1,  which  was  used  to  estimate  the  required  MTTF  or 
scale  parameter,  6,  of  the  failure  generating  function  for  this  simula¬ 
tion.  Once  the  shape  and  scale  parameters  have  been  estimated  for  a 
given  equipment/component  failure  distribution,  then  the  MTTF  can  be 
estimated  by  using  the  table.  For  example 


=  6  •  r  (i  +  i/K) 

Mean  Failures  Per  Hr  = _ i _  = _ \ _  = _ V 

E(t)  0*r(l+l/K)  MTTF 

From  the  above,  given  the  MFPH  and  K,  the  scale  parameter 
can  be  estimated  from: 

_ 1 _  MTTF 

9  =  MFPH  •  r(l  +  l/K)  =  r  (1+1/K) 

Thus,  by  estimating  the  shape  parameter,  K,  and  determining 
one  of  the  other  two  parameters  Q  or  MFPH,  we  can  derive  an  estimate 
for  the  remaining  variable  when  the  failure  distribution  is  of  the  Weibull 
form: 

f(t)  =  tk_1  e 
0k 

The  table  is  constructed  with  9  =  1  and  the  K  parameter  is 
varied  from  0  to  10  in  increments  of  .  25.  All  that  is  necessary  to 


determine  the  MFPH  cr  Q  value  is  to  obtain  the  tabulated  mean  value. 
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TABLE  5 

STANDARDIZED  TABLE  TO  DETERMINE  THE  STEADY 
STATE  MFPH  OF  COMPONENTS  WHICH  HAVE  A 
WEIBULL  TIME  TO  FAILURE  DISTRIBUTION 


SHAPE 


STANDARDIZED  TABLE 

EXPECTED  MEAN 


0.250 

0.04166667 

0.500 

0. 50000000 

0.  750 

0. 83988503 

1.000 

1. 00000000 

1.250 

1. 07367142 

1.  500 

1.  10773224 

1.  750 

1.  12281597 

2.000 

1. 12837909 

2.  250 

1.  12901184 

2.  500 

1.  12706082 

2.  750 

1.  12377663 

3.  000 

1. 11984669 

3.  250 

1.  11565629 

3.  500 

1. 11142324 

3.  750 

1. 10726975 

4.000 

1.  10326244 

4.  250 

1. 09943514 

4.  500 

1.  09880223 

4.  750 

1. 09236656 

5.  000 

1. 08912425 

5.  250 

1. 08606765 

5.  500 

1. 08318715 

5.  750 

1. 08047223 

6.  000 

1. 07791221 

6.  250 

1. 07549660 

6.  500 

1.  07321534 

6.  750 

1.  07105894 

7.000 

1. 06901857 

7.  250 

1.  06708595 

7.  500 

1.  06525347 

7.  750 

1. 06351410 

8.  000 

1.  06186141 

8.  250 

1. 06028937 

8.  500 

1.  05879261 

8.  750 

1. 05736604 

9.000 

1. 0560051 1 

9.  250 

1.  05470550 

'  .  500 

' 1. 05346334 
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TABLE  5  (Continued) 


SHAPE 


EXPECTED  MEAN 


9.  750  1.05227505 

10.000  1.05113725 


given  the  estimated  K  parameter  and  then  divide  the  value  by  the  pre¬ 
determined  scale  parameter,  or  MFPH  respectively. 

Examples: 

(1)  Given  K  =  3.  5 
and,  t1  =  950. 

Obtaining  from  the  table,  MEAN  =  1.11142324 

the  MFPH  =  MgftN  =.0011699 
950 

(2)  Given  K  =  2 

and,  MTTF  =  1500 
or,  MFPH  =  (1 500)" 1 

Obtaining  from  the  t.  ble,  MEAN  =  1.  12837909 

e=  =1692.5685 

(1500f  1 


In  the  case  of  the  consta.it  failure  rate,  the  shape  parameter, 
k,  for  the  Weibull  cumulative  density  function  is  1  and  the  scale 
parameter,  9  ,  is  directly  proportional  to  the  MTBF  as  follows: 


r  * 

MTTF  =  1  t  f(t)  dt 

J0 


=/*  fa  'k'‘  e’l»j  ',t'  K>0’e>0 


e 
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=  dr  (i  +  i/K) 

With  K  =  1,  MTTF  becomes: 

MTTF  =  Or  (1  +  l/l) 
or 

MTTF  =  Or  (2) 

But  r( 2)  is  equal  to  1,  therefore 

MTTF  =  6 *  1  =  6 

and, 

MFPH  = - l- -  =  I 

MTBF  6 

Thus,  when  the  underlying  failure  distribution  follows  the 
exponential  distribution,  we  can  equate  the  scale  parameter  to  the 
reciprocal  of  the  MFPH,  or  equate  the  scale  parameter  directly  to  the 
MTTF.  Under  these  conditions,  we  can  predict  the  steady  state  MFPH 
by  computing  the  reciprocal  of  the  MTTF  or  of  8. 

Summary 

Of  the  three  cases  investigated,  only  the  constant  failure  rate 
case  appeared  to  be  in  steady  state  immediately.  The  decreasing  (DFR) 
and  increasing  (IFR)  failure  rate  cases  result  in  significant  deviations 
from  the  failure  rates  predicted  by  the  present  Air  Force  methods 
which  may  lead  to  supply  difficulties  in  supporting  weapons  systems 
during  the  transition  period. 


The  transition  period  prior  to  the  attainment  of  steady  state 
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conditions  for  the  DFR  and  IFR  exceeds  12  months  for  the  MTBF's  and 
parameters  tested  and  may  actually  be  longer  than  24  months  for  large 
MTBF's.  Systems  in  transition  tend  to  approach  steady  state  more 
rapidly  for  lower  MTBF  values,  and  the  random  relative  fluctuations 
appear  to  be  less  erratic  than  those  for  larger  MTBF  values. 

The  plots  of  the  system  failure  rates  indicate  that  upon  attain¬ 
ment  of  steady  state,  the  expected  future  failure  rates  for  the  system 
can  be  approximated  or  predicted  by  the  present  Air  Force  methods 
quite  satisfactorily.  In  the  limit,  the  population  failure  distributions 
appear  to  exhibit  the  random  oscillations  around  the  steady  state  fail¬ 
ure  rate  characteristic  of  the  constant  failure  rate  case. 


CHAPTER  VII 


CONCLUSION 


Our  original  hypothesis  that  the  fleetwide  distributions  of  fail¬ 
ures  for  any  replaceable  equipment/ component  is  at  all  times  a  Poisson 
process  and  that  the  assumption  of  a  constant  failure  rate  is  always 
valid  in  the  estimation  of  future  expected  failures  is  not  supported  by 
the  simulation  results  and  we  reject  the  hypothesis. 

The  simulation  results  using  the  Weibull  failure  generating 
function  as  the  underlying  failure  distribution  for  the  population  of 
equipments/ components  indicate  that  a  deviation  between  the  simulated 
failures  and  those  failures  predicted  by  a  Poisson  process  is  to  be 
expected  when  aircraft  components  are  introduced  incrementally  into 
the  Air  Force  inventory.  The  transition  period  between  the  initial  level 
of  failures  and  the  steady  state  level  of  failures  is  related  to  the  MTTF 
of  the  introduced  components  and  the  characteristic  shape  parameter 
of  the  underlying  failure  distribution. 

Since  the  shape  parameter  of  the  Weibull  failure  distribution 
determines  the  general  "shape"  of  the  distribution,  we  can  approximate 
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the  exponential  distribution  and  the  normal  distributions  which  are  the 
basic  statistical  forms  assumed  to  describe  the  failure  characteristics 
of  components.  As  the  shape  parameter  is  increased  above  the  value  1, 
the  Weibull  tends  to  approach  a  normal  distribution.  The  larger  the  K 
parameter  becomes  the  smaller  the  deviation  from  the  predicted  failure 
rate  becomes,  until  in  the  limit,  i.  e.  ,  oo  ,  the  distribution  has  no 
deviation  and  the  failures,  MFPH,  can  be  approximated  by  ' 

This  can  be  seen  from  the  following: 

1  1 

MFPH  =  MTTF  =  dr  (1  +  l/k) 

lim  r(l  +  1/K)  =  r(l)  =  1 

K-*oo 

therefore, 

lim  MFPH  =  lim  1  =  i  =  1 

K->oo  K->  oo  Qr  (1+1  /K)  6  MTTF 

As  can  be  seen  from  the  Standardized  Table  (Table  5)  ar.d  graph 
(Graph  1)  the  maximum  average  deviation  from  the  predicted  exponential 
MFPH  the  range  2.  <  K  <  2.  25.  The  maximum  deviation  occurs  at 
approximately  K  =  2.  17.  Such  a  K  value  is  quite  possible  in  actual 
experience  with  equipment/component  failure  rates,  since  such  a  K 
value  places  the  underlying  failure  distribution  "between"  the  exponential, 
i.  e.  ,  K  =  1,  and  an  approximation  to  the  shape  of  the  normal  distribu¬ 
tion,  i.  e.  ,  K  =  3.  Thus,  when  the  Air  Force  acquires  equipments/ 
components  incrementally  and  the  actual  underlying  failure  distribution 
is  of  the  increasing  failure  rate  type,  we  should  not  be  surprised  to 
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observe  increasing  failure  rates  during  the  transition  period  to  steady 
state. 

When  we  have  equipments /components  which  are  characterized 
by  decreasing  failure  rates,  the  Weibull  failure  generating  function  has 
a  shape  parameter  0  <  K<  1.  Such  a  case  has  been  explored  in  this 
simulation  and  the  results  have  been  indicated  in  the  Analysis  Section 
of  this  report  as  the  Decreasing  Failure  Rate  Case.  Simulation  results 
indicate  that  a  shape  parameter  value  of  .  5  causes  the  experienced  fail¬ 
ure  rate  to  approach  the  predicted  failure  rate  asymptotically  from 
above.  The  deviation  from  the  predicted  failure  rate  and  the  simulated 
failure  rate  increasing  directly  with  the  MTBF.  In  such  a  situation,  we 
should  expect  to  see  larger  initial  failure  rates  which  decrease  through¬ 
out  the  introduction  period  and  eventually  settling  down  when  steady  state 
is  reached. 

In  the  case  of  constant  failure  rates,  i.  e.  ,  a  shape  parameter 
K  =  1,  the  simulation  results  indicate  that  failures  will  be  randomly 
distributed  about  the  predicted  failure  rate  increasing  in  absolute 
deviation  from  the  predicted  failure  rate  as  the  MTBF  increases.  We 
should  not  be  surprised  nor  concerned  with  apparent  increases  or 
decreases  in  the  number  of  failures  observed  when  the  components  have 
an  underlying  exponential  failure  distribution,  since  such  deviations 
from  the  expected  or  predicted  failure  rate  are  due  to  random  fluctu¬ 
ations. 


The  simulation  results  also  indicate  that  after  steady  state  has 
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been  attained  for  a  system  of  equipments /components  the  observed 
distribution  of  failures  resembles  the  random  fluctuations  of  the 
constant  failure  rate  case.  This  is  the  result  proved  by  Khintchine  as 
previously  indicated  in  Chapter  III.  Thus,  it  appears  that  upon  the 
attainment  of  steady  state,  the  expected  failures  can  be  estimated  by 
assuming  a  Poisson  process,  regardless  of  the  underlying  failure  gene¬ 
rating  distribution;  and,  in  steady  state  the  MTBF  should  equal  the 


MTTF. 


APPENDIX  A 


COMPUTER  SIMULATION  PROGRAM 


The  program  developed  for  this  research  is  shown  in  Figure  4 
and  a  simplified  logic  chart  is  shown  in  Figure  5.  This  program  is  a 
FORTRAN  language  program  using  the  GASP2B  simulation  package 
developed  by  A.  Alan  B.  Pritsker  and  Philip  J.  Kiviat.  (16)  The  use 
of  GASP2B  simulation  package  greatly  simplifies  the  programming 
required  to  be  written  for  this  study. 

Basically  the  program  is  composed  of  a  Main  Program  and  six 
Subroutines  as  shown  in  Figure  4.  The  Main  Program  is  used  to 
initialize  non-GASP  variables  and  to  call  the  GASP2B  program  package. 
When  the  GASP  package  is  called,  the  internal  GASP  executive  routines 
control  all  further  simulation  everts  until  simulation  termination  is 
reached. 

The  GASP  executive  routines  call  Subroutine  EVNTS  which  in 
turn  selects  the  appropriate  INTRO,  FAILR,  REPORT,  or  ENDSM 
subroutines  based  upon  the  currently  processed  event. 

The  INTRO  subroutine  is  primarily  used  to  introduce  the 


Preceding  gage  blank 
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aircraft  into  the  inventory  incrementally  until  the  total  fleet  size  has 
been  attained.  Each  unit  introduced  incrementally  to  the  system  is 
assigned  a  time  to  failure  (service  time)  and  internally  filed  in  a  future 
events  file. 

The  FAILR  subroutine  is  called  by  the  GASP  routine  whenever  a 
failure  event  occurs.  To  simulate  the  instantaneous  repair/replacement 
of  failed  components,  this  subroutine  immediately  computes  a  new  time 
to  failure  and  then  refiles  this  new  time  to  failure  in  the  future  events 
file.  Statistics  are  collected  on  the  service  time  of  each  failure. 

The  REPORT  subroutine  schedules  the  next  monthly  reporting, 
and  records  the  total  period  failures  and  service  times  which  have  been 
experienced.  A  check  is  made  against  the  total  simulation  time  for 
each  run  and  the  simulation  run  is  terminated  if  required. 

The  ENDSM  subroutine  sets  specific  GASP  variables  for  the 
termination  of  the  current  simulation  run.  A  check  is  made  to  deter¬ 
mine  if  the  total  number  of  simulation  runs  has  been  completed  and  the 
OUTPUT  subroutine  is  called  as  required. 

Tne  OUTPUT  subroutine  has  been  written  to  provide  the  total 
simulation  output  when  the  total  number  of  simulation  runs  has  been 
completed.  In  this  research,  a  simulation  run  is  approximately  4  years 
long  and  each  simulation  run  is  performed  40  times.  The  OUTPUT 
subroutine  then  takes  an  average  of  all  statistics  collected  over  the  40 
runs  before  the  output  is  printed.  In  this  manner,  we  minimize  the 

effect  of  extreme  values  in  our  output  data.  The  printed  output  format 
and  data  arc  shown  in  the  simulation  output  in  Appendix  B. 


Figure  4.  --Computer  Program 


Figure  4  continued 
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APPENDIX  B 


SIMULATION  OUTPUT  DATA 


This  appendix  contains  the  final  simulation  output  listings  and 
plots  for  the  various  parameters  investigated  in  this  study.  Each  list¬ 
ing  contains  the  average  number  of  failures  per  period,  the  period 
operating  hours,  the  average  number  of  failures  per  operating  hour  per 
period  and  the  period  standard  deviation  of  failures  per  operating  hour. 
The  legend  at  the  bottom  of  each  listing  describes  the  conditions  for 
which  the  data  was  generated.  Definitions  of  these  conditions  are  as 
follows : 

NUMBER  IN  FLEET- -the  maximum  size  of  the  fleet  to  be  in 
the  inventory. 

INTRO  INTER  VAL- -the  time  interval  between  increments  of  the 
fleet  in  hours.  (120  hrs.  equals  1  Report  Interval). 

REPORT  INTERVAL- -the  time  length  of  each  reporting  period 
in  hours. 

SCALE-  -a  parameter  of  the  Weibull  distribution  determined  by 
the  equation 
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MTBF 


SHAPE- -the  other  parameter,  K,  of  the  Weibull  distribution 
tested  in  this  study. 

MTBF- -the  Mean  Time  Between  Failures  under  the  assumption 
of  an  exponential  time  to  failure  distribution. 

NUMBER  INTRO/PERIOD- -the  number  of  items  introduced  to 
the  simulation  until  the  total  fleet  size  is  attained. 

The  failures  per  ops  hour  have  been  plotted  for  each  listing  and  are 


included  in  Chapter  VI  of  this  report. 

The  listings  are  presented  as  follows: 


TABLE  NO. 
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SHAPE 

B-l 

500 

.  5 

B-2 

1000 

.5 
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.  5 

B-4 

50 

1.  0 
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100 

1.0 
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500 

1.  0 
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1.  0 
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1.0 
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50 

2.  0 

B-l  0 

50 

3.  0 

B-l  1 

50 

4.  0 
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5.  0 
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2.  0 

B-14 

100 

3.  0 

B-15 

100 

4.  0 

B-16 

100 

5.  0 
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1.  5 
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2.  0 

B-19 
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3.  0 

B-20 
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4.0 

B-2 1 

500 

5.  0 

B-22 

1000 

1.  5 

B-23 

1000 

2.  0 
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Figure  15 
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Figure  32 
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